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Abstract 

A Laplacian matrix, L — (iij) £ R"^", has nonpositive off-diagonal entries and zero row sums. 
As a matrix associated with a weighted directed graph, it generalizes the Laplacian matrix of an 
ordinary graph. A standardized Laplacian matrix is a Laplacian matrix with — ^ < ^ij < at 
j 7^ i. We study the spectra of Laplacian matrices and relations between Laplacian matrices and 
stochastic matrices. We prove that the standardized Laplacian matrices L are semiconvergent. 
The multiplicities of and 1 as the eigenvalues of L are equal to the in-forest dimension of the 
corresponding digraph and one less than the in-forest dimension of the complementary digraph, 
respectively. We localize the spectra of the standardized Laplacian matrices of order n and study 
the asymptotic properties of the corresponding domain. One corollary is that the maximum 
possible imaginary part of an eigenvalue of L converges to ^ as n oo. 

1 Introduction 

We consider the matrices L = {£ij) E K."^" such that 

4j < 0, j^i, (1) 

n 

= 0, i = l,...,n. (2) 

If, in addition, L is symmetric and £ij G {0,-1}, j ^ i, then L is the Laplacian matrix of an 
undirected graph (see, e.g., ^): £ij = —1 iff the graph has an edge between vertex i and vertex 
j, j ^ i; £ii — the degree of vertex i. For the classical results on Laplacian matrices, we refer to 

o nsi nni n cn in ini CHI d ini . 

In the same way, every matrix of the form of ^ and jSJl can be associated with a weighted 
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directed graph T whose arc weights are strictly positive: 

minus the weight of arc \i j ^ i and F has an (i, j) arc; 

t^j ^ I 0, if j 7^ i and r has no (i, j) arc; (3) 

(cf. the KirchhofF matrix 

According to the famous matrix-tree theorem for digraphs the cofactor of such a 

matrix i(r) is equal to the total weight of spanning trees converging to i in F, where the weight 
of a subgraph is the product of the weights of its arcs. Generalizations of this theorem have been 
obtained in jJUj and [SJ |3] . We call the matrices that satisfy and (0) Laplacian matrices. 

It is easily seen |3J p. 258] that every Laplacian matrix is a singular M-matrix (for a detailed 
discussion of M-matrices we refer to ,2,). By Gersgorin's theorem (see, e.g., HSj), the real part of 
each nonzero eigenvalue of L is strictly positive. 

Laplacian matrices are closely connected with stochastic matrices. The basic relation is that the 
Laplacian matrices can be defined as a(/ — P) matrices, where a > 0, / is the identity matrix, and P 
is a stochastic matrix. That is why the Laplacian matrices frequently appear in the theory of Markov 
chains.^ Dmitriev and Dynkin \7\ note that the matrices — L, where L are Laplacian matrices, are 
the matrices of transition probability densities of continuous time Markov chains. They strengthen 
the aforementioned result on the nonnegativity of the real parts of the Laplacian eigenvalues by 
proving that for every eigenvalue A, 

In this paper, we further characterize the spectra of the Laplacian matrices. To achieve a better 
comparability of the results for various n, we consider standardized Laplacian matrices. These are 
the Laplacian matrices L = {£ij ) with 

--<i^J<0, J ^i, i = l,...,n. (5) 
n 

For such matrices, 

0<£u<l--, i = l,...,n. (6) 

n 

If the class Qb of weighted digraphs with positive arc weights not exceeding 6 > is considered 
and L(F) is the Laplacian matrix of a weighted digraph F on n vertices, then the standardized 
Laplacian matrix associated with F in this class is, by definition, i(F) = (nb)^^ L{T). 

It follows from 1, Theorem 1] (or from the earlier results (2.7) and (2.10) of JS]) that, in the 
case of symmetric Laplacian matrices, the eigenvalues of L belong to [0, 1]. In the general case, (0}, 
and Gersgorin's theorem immediately provide 

^Studying the spectra of nonsymmetric Laplacian matrices is also needed for solving control problems that involve 
graphs 0. 
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Proposition 1 The spectrum of a standardized Laplacian matrix belongs to the meet of the closed 
disk with radius l — l/n centered at 1 — l/n and the closed smaller angle bounded by the two half-lines 
drawn from through the points e~'^^~^^^ and e*-^~"'*. 

For n — 7, Proposition ^ is illustrated by Figure 1. 





.Im(A) 

n //.__..__..__.__.. \ 






Re(A) 









V n — 1 \ 
\ n 













Figure 1. 



In Section 2, we show that the eigenvalues of L are also the eigenvalues of a certain stochastic 
matrix and that L is semiconvergent; the multiplicities of and 1 as eigenvalues of L are related 
with the in- forest dimensions of the corresponding digraph F and the complementary digraph Fc, 
respectively. In the subsequent sections, we advance in localizing the spectra of the nonsymmetric 
Laplacian matrices. 



2 Laplacian matrices and stochastic matrices 

Let 

K = I-J, (7) 

where J £ M"^" is the matrix with all entries 1/n. Obviously, K E R"^" is the standardized 
Laplacian matrix of the complete digraph with all arc weights b in the class Gf 
Define the matrices 

P = L + J (8) 

and 

Lc^K-L. (9) 
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P is stochastic, since it is nonnegative with all row sums 1. In the class Qb, Lc is the standardized 
Laplacian matrix of the complementary weighted digraph where the weight attached to arc (i, j), 
j 7^ I, is 6 — Wij, Wij being the weight of this arc in F. If Wij = b then Fc has no (i,j) arc and vice 
versa: if F has no arc, then the weight of in Fc is b. It follows from Q, and ^ that 

P = I-L,. (10) 

Let <t{A) be the spectrmii of a square matrix A. 

Theorem 1 Suppose that L is a standardized Laplacian matrix, P is the stochastic matrix defined 
by @, and Lc is the matrix defined by @l. Then for A ^ {0,1}, the following statements are 
equivalent: 

(a) A e a-(Z), 

(b) A e cr(P), and 

(c) (1 - A) e a(Z,), 

and these eigenvalues have the same geometric multiplicity. Furthermore, v is an eigenvector of L 
corresponding to \ ^ {0, 1} if and only if the vector x such that 

I-T^^^ (11) 



1 - ^, 

is an eigenvector of P corresponding to A and of Lc corresponding to \ — X. 

Proof. Let v be an eigenvector of L corresponding to an eigenvalue A ^ {0, 1}. Since LJ = 
and = J, we have 

Conversely, let x be an eigenvector of P corresponding to A ^ {0, 1}. Observe that (|ll|l is 
equivalent to 

v^(l-Ox. (12) 



A 



whenever A ^ {0, 1}. Indeed, since ~ J, 



I-aJ^ x) ^ X 1-A^A(1-A)~' 
Using LJ = 0, we get 

Lv = L — — ^ X = Lx = [P — J)x = Xx — Jx — x(^ — — ^ x = Xv. 

The assertions about Lc follow from (|10|l . The geometric multiplicities of the eigenvalues are the 
same, since the transformation Hll|l is nondegenerate: p2|l defines its inverse. □ 
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Theorem 2 Let fj^{X), fp{X), and fj^ (A) be the characteristic polynomials of L, P, and Lc, respec 
tively. Then for all X ^ {0,1}, 

X - 1 



/p(A) 



X 

i-iy 



A 



1 - A 



/i(l-A). 



Proof. Let c^^\ . . . , c*^"^ be the columns of the characteristic matrix XI — L. Let 1 = {-, . 
Then, by ®, 



(13) 
(14) 

2_\T 
n / ■ 



/p(A) = det 
= det 
- det 



(2) 



. ,c 



det 



det 



(15) 



Note that det [c^'^',c' 



) ^(2) 



/r(A) and 1 = (An)-i(c(i) + • • • + €(")), provided that A 7^ 0. 



Therefore 



det 



l,c' 



(2) 



.,c' 



(") 



det 



An 



fe=2 



Ai 



This expression works for each of the last n items of (|15f) . Substituting this in (|15|) gives /p(A) = 
/^(A) -n3i;;/£(A) = ^/£(A). To prove (0, note that by COI), for aU A 7^ 1, /^^(A) =det(A/ - / + 
P) = (-1)" det(/ -XI-P) = (-l)"/p(l - A) = (-!)"-! ^/~(1 - A). □ 

Remark 1. Since H13() and H14() are true for all A, except for and 1, they become valid for all A 
after reducing the common factors A and 1 — A, respectively. 

An analogy of (|14|l for symmetric matrices is due to Kelmans Property 5] . 
Recall that a square matrix A is semiconvergent if limfe^oo A^ exists. The index of A, indA, is 
the least fc = 0, 1, . . . , rt — 1 such that rankv4'^+^ — rankyl'^. 

Theorem 3 For every standardized Laplacian matrix L and the corresponding stochastic matrix P 
defined by L and P are semiconvergent. 

Proof. Let p{A) be the spectral radius of a matrix A. A square matrix A is semiconvergent 
(see, e.g., |21 Problem 4.9]) if and only if it satisfies: 

(a) piA) < 1, 

(b) If A e a{A) and |A| = 1, then A = 1, and 

(c) rank(/ - A)^ = rank(/ - A), i.e., ind(/ - A) e {0, 1}. 

We now prove that P satisfies (a), (b), and (c). (a) results from the stochasticity of P. 

(b) Since all diagonal entries of the singular M-matrix Lc — I — P are less than 1, |A| — 1 implies 
A = 1 (see, e.g., page 153]). 

(c) Since the index of any Laplacian matrix is 1 0] Proposition 12], we have rank(/ — P)^ = 
rank(/ — P), and P satisfies (c). Thus, P is semiconvergent. 
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Finally, it follows from © and LJ^O that L = P - J , = L{P - J) = LP = {P - J)P, 
^ L{P - J)P'=-2 ^ £pfc-i = (p - J)P''-^, .... 

Therefore, 

lim ^{P- J) lim P''"\ (16) 

— ^oo k — *oo 

and L is semiconvergent as well. □ 

Let 771^ (A) be the algebraic multiplicity of A as an eigenvalue of A\ let Va{X) be the set of 
eigenvectors of A corresponding to A. An eigenvalue is called semisimple if its algebraic and geometric 
multiplicities coincide. The in-forest dimension of a digraph is the minimum possible number of 
trees in a spanning converging forest (also called in-forest) of the digraph ^ p. 255]. 

Theorem 4 Let d and dc be the in-forest dimensions of the digraph corresponding to L and the 
complementary digraph, respectively. Then: 

(i) ?7i2^(0) = d, mj^(l) — dc — 1; 

(ii) TOp(O) = d - 1, mp(l) = dc] 

(iii) (1) = — 1, m£ (0) = dc, 
and these eigenvalues are semisimple. 

(iv) Ifve ViiO) and Kv ^ 0, then Kv € Vp{0) = Vi^{l); 
ifxe Vp{l) = VijQ) and Kx ^ 0, then Kx G Vi{l). 

Proof. By 0] Proposition 12], rankL — n — d and indi = 1, therefore m^(0) = d and is a 
semisimple eigenvalue. By Theorem |3| L is semiconvergent. Thus, L satisfies condition (c) in the 
proof of Theorem 13 which implies the semisimplicity of 1, provided that 1 G cr{L). Indeed, by (c), 
the maximum order of the Jordan blocks corresponding to 1 in the Jordan form of L is at most 1. 
By virtue of 1)14(1 and Remark 1, rn-j^{l) = mj^ (0) — 1. The statements regarding Lc follow from 
those about L; after that, (ii) and the assertions regarding P result from (|10|l . The proof of (iv) is 
straightforward. □ 

For the case of symmetric Laplacian matrices, the identities mj^{Q) — d and 771^^(1) = dc — 1 
are due to Kelmans |14 |. In that case, d and dc reduce to the numbers of components of the 
corresponding graphs. 

In the following sections, we focus on the eigenvalues with nonzero imaginary parts. 



3 A region that contains the Laplacian spectra 

Theorem 5 Let L be a standardized Laplacian matrix of order n. Then the spectrum of L belongs 
to the meet of: 

- two closed disks, one centered at l/n, the other centered at l — l/n, each having radius 1 — 1/n, 

- two closed smaller angles, one bounded with the two half-lines drawn from 1 through g^^Tri/n 
and e^'^*/", the other bounded with the half-lines drawn from through e"^?"^")' and 6^^^"-*% and 

- the band I Im(z) I < 77- cot 7?- ■ 
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Proof. In addition to Proposition^ note that the domain specified in that proposition contains 
the spectrum of too. Therefore, the spectrum of P — I ~ and, by Theorem ^ the spectrum 
of L belongs to the disk centered at 1/n and the angle with vertex at 1, as specified in Theorem [S] 

Furthermore, according to Pick's theorem PHj, if A is an eigenvalue of A € R"^", then | Im(A)| < 
g cot(7r/27i), where g — max|cfej| and C — (ckj) ~ ^{A ~ A*). If A — L, then 5 < ^ by virtue 
of ISl). This completes the proof. □ 

Remark 2. It follows from the stochasticity oi I — L and [S] Theorem 4] that a{L) contains an 
eigenvalue with argument f — ^ if and only if the corresponding digraph is a Hamiltonian cycle on 
n vertices. In this case, such an eigenvalue A is unique, |A| < ^sin^, and Im(A) < — sin — . These 
bounds are reached when L — ^{I ~ Q), where Q is the standard circulant permutation matrix 
(see the following section) . The components of any corresponding eigenvector are the vertices of a 
regular polygon. Similarly, a{L) contains an eigenvalue that belongs to [Ije'^] if and only if the 
complementary digraph is a Hamiltonian cycle on n vertices. As above, such an eigenvalue A' is 
unique and Im(A') < ^ sin ^ZL. 

TheoremEland Remark 2 are illustrated by Figure 2, where n — 7. 




Figure 2. 

The region R specified by Theorem |31 (hatched in Figure 2) is a hexagon whenever 4 < n < 18, 
however its boundary contains arcs when n > 18. This region reduces to the [0, 1] segment when 
n = 2 and to the rhombus with vertices 0, (0.5, (2'\/3)^^), 1, and (0.5, —{2^)^^) when n = 3. 
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4 The polygon of Laplacian eigenvalues 



In this section, we introduce a polygon that consists of Laplacian eigenvalues. As above, the order 
n of standardized Laplacian matrices is fixed. 

Let Q — (qkj) be the standard nxn circulant permutation matrix: qkj = 1 if .7 — fc G {1, 1 — n} 
and qkj — otherwise. Consider n — 1 matrices 

Lk^-(kI-Q-Q^ Q'^), k = l,...,n-l. (17) 

Since Q*, z = 1, . . . , n— 1, are permutation matrices whose patterns of ones are different entrywise, 
the Lfc's have zero row sums, and their off-diagonal entries belong to [— i,0]. Thus, they are 
standardized Laplacian matrices. 

The characteristic polynomial of Q is /q(/^) = ^" — 1, therefore the spectrum of Q contains 

[1 = e-2W". (18) 

Consequently, the spectrum of contains 

Afc = ^fe) (19) 

n 



i (/c - g-^'r*/" g-2fc^,:/„j ^ (20) 



as well as A^. 

By means of the identity 



sin ^ 

^ai _|_ g2ai ^ y ^kai ^ 2 g(fc+l)Qj/2 



(I20II simplifies to: 

Xk^n-^(k-'^ e-(-^)-/"l ^ A. _ ^ (^^^ i±±l> _ ^ ^.^ (fc±l)^ 

n n sin ■ 



(21) 



Let S be the convex polygon with 2{n — 1) vertices: Aq = 0, Ai, . . . , A„_2, A„_i = 1, A„_2, . . . , Ai. 

Theorem 6 Every point of the polygon S is an eigenvalue of some standardized Laplacian matrix 
of order n. 

Proof. As shown above, every vertex A^ of S is the eigenvalue of the corresponding Lk- Every 
other point s G 5* (in what follows, we suppose Im(s) > 0, since the spectrum of a real matrix is 
symmetric with respect to the x-axis) can be represented by a convex combination /(O, Afe,Afc+i) 
of three vertices of S: and the two neighboring vertices Afc and Afc+i such that the half- line 
[0, s) meets the segment [Afc,Afc-|-i]. By definition, Afc and Afc+i are the values of the polynomials 
(|19|l for /i (we denote these polynomials by p(-) and q(-), respectively); furthermore, they are the 
eigenvalues of the matrices Lk — p{Q) and Lk+i — q{Q), respectively. Thus, s = f{0,p{^j,),q{^)), 



8 



where /(0,p(-), is a polynomial. Consequently, s is the eigenvalue of the matrix f{0,p{Q), q{Q)), 
which is a standardized Laplacian matrix as a convex combination of such matrices. □ 

It is easily seen that the polygon S reduces to [0, 1] when n = 2. For n = 3, S is the rhombus 
with vertices 0, (0.5, (2\/3)^^), 1, (0.5, — (2\/3)^^), and it coincides with the region described in 
Theorem |S1 as well. 

Let 

z{n) = sup{| Im(A)| : A is an eigenvalue of a standardized nxn Laplacian matrix}. (22) 
By TheoremEl for every n = 2, 3, . . . , z{n) < ^ cot 
Proposition 2 If n is odd, then 

moreover, z{n) = Im(A(„_i)/2): where A(„_i)/2 *s defined by (jj.^jl . 
The proposition is proved by substituting fc = (n — l)/2 in (|21(l . 

If ri > 2 is even, then the upper bound cot ^ is not reached by the imaginary parts of the 
vertices of S; more specifically, the following proposition is easy to verify. 

Proposition 3 If n > 2 is even, then 

max Im(Afe) = Im(A„/2) = — cot — < — cot — . 

The following conjecture is yet unproved. 

Conjecture 1 The set of all eigenvalues of standardized Laplacian matrices of order n coincides 
with the polygon S . 



5 Some asymptotic properties of the Laplacian spectra 

Theorem 7 The boundary of the polygon S converges, as n ^ oo, to the curve made up by the 
parts of two cycloids whose parametric equations are: z{t) — x{t) +iy(T) and z{t) — x{t) ^iy^r), 
where 

x{t) = (27r)-i(T- sinr), (24) 
2/(t) = (2^)-i(l-cosr), (25) 

and T £ [0, 27r] . 

Proof. For each r G [0, 27r), consider any positive integer sequence k{n) such that lim -^^^ = 
^. Then, by virtue of lEH), 

lim Afc(„) = — (t - 2 sin - cos - + 2isin^ - ) = — (t - sinr + i(l - cost)). 
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Taking the real and imaginary parts of this expression, we obtain the desired statement. □ 

Figure 3 shows the polygons S" at n = 4 and n — 5, along with the limit curve specified in 
Theorem 
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Figure 3. 



Corollary 1 lim z{n) = I/tt, and the convergence is from below. 

n — *oo 

This corollary is easily proved using PropositionsElandOl 
Acknowledgments. We are grateful to two anonymous referees for their helpful comments. 
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